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Triacylglycerol (TAG) accumulates in plant seeds as a major renewable source of carbon for 
food, fuel and industrial feedstock. Approaches to enhance TAG content by altering lipid 
pathways and genes in vegetative parts have gained significant attention for biofuel and 
other applications. However, consequences of these modifications are not always studied 
in detail. In an attempt to increase TAG levels in leaves we previously demonstrated 
that a novel substrate, monoacylglycerol (MAG), can be used for the biosynthesis of 
diacylglycerol (DAG) and TAG. Transient expression of the Mus musculus monoacylglycerol 
acyltransferases MGATl and 2 in the model plant Nicotiana benthamiana increased 
TAG levels at 5 days post-infiltration (dpi). Here we show that increased TAG and 
DAG levels can be achieved as early as 2 dpi. In addition, the IVIGATl infiltrated 
areas showed senescence-like symptoms from 3 dpi onwards. To unravel underlying 
molecular mechanisms, lllumina deep sequencing was carried out (a) for de-novo 
assembling and annotation of N. bentliamiana leaf transcripts and (b) to characterize 
A//G/\r?-responsive transcriptome. We found that /WG/Ar7-responsive genes are involved 
in several processes including TAG biosynthesis, photosynthesis, cell-wall, cutin, suberin, 
wax and mucilage biosynthesis, lipid and hormone metabolism. Comparative analysis 
with transcript profiles from other senescence studies identified characteristic gene 
expression changes involved in senescence induction. We confirmed that increased TAG 
and observed senescence-symptoms are due to the MAG depletion caused by MGATl 
activity and suggest a mechanism for IVIGATl induced TAG increase and senescence-like 
symptoms. The data generated will serve as a valuable resource for oil and senescence 
related studies and for future N. benthamiana transcriptome studies. 
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INTRODUCTION 

Demand for production and use of biofuels is rapidly increasing 
and current agricultural practices are unable to produce sufficient 
amounts cost-effectively. Multiple approaches are being employed 
to develop plants with altered carbon and energy content to meet 
this demand. Accumulation of energy-dense lipids like triacyl- 
glycerols (TAG) in the vegetative parts of plants has attracted sig- 
nificant attention due to the potential to produce large amounts 
of energy rich biomass (Vanhercke et al., 2014a,b). Approaches to 
increase TAG levels in diverse plant tissues have targeted differ- 
ent steps of TAG biosynthesis and storage that can be generally 
divided into (a) "Push"- mobilizing carbon flux to increase fatty 
acid biosynthesis, (b) "Pull" — TAG assembly, and (c) "Protect" — 
increasing TAG storage/decreasing TAG breakdown (Vanhercke 
et al, 2013a,b). Examples of the "Push approach" include ectopic 
expression of transcription factors (LECl, LEC2 and WRIl) that 



control seed development and maturation (Santos Mendoza et al., 
2005; Andrianov et al., 2010; Liu et al., 2010) as well as diverting 
carbon flow from starch to lipid biosynthesis by silencing ADP- 
glucose-pyrophosphorylase (AGPase) in plant leaves (Sanjaya 
et al, 2011 ). "Pull approaches" consist of overexpression of acetyl- 
CoA carboxylase (Klaus et al., 2004) and overexpression of dia- 
cylglycerol acyltransferase (DGAT) and other enzymes involved 
in TAG assembly (Bouvier-Nave et al, 2000; Jako et al, 2001; 
Weselake et al., 2008; Petrie et al, 2012). Finally, retaining fatty 
acid (FA) content ("Protect approach") by blocking FA move- 
ment to peroxisomes or via other mechanisms can also increase 
TAG levels in leaves (Xu et al, 2005; Slocombe et al, 2009). Many 
studies which combine several of these approaches (Andrianov 
et al, 2010; Sanjaya et al, 2011; Vanhercke et al, 2013a) have 
provided "proof of concept" results for oil increase, although the 
TAG levels achieved are not yet necessarily of industrial relevance. 



www.frontiersin.org 



May 2014 | Volume 5 | Article 204 | 1 



Divi et al. 



Monoacylglycerol-mediated triacylglycerol synthesis and senescence 



However, recently substantial progress had been made in this 
direction by metabolic engineering of plant leaves for high energy 
dense oil-seed like TAG (Vanhercke et al., 2014a). 

In plants diacylglycerol (DAG), the immediate precursor 
of TAG, is synthesized mainly by: (1) Kennedy pathway 
for d& novo DAG synthesis (Kennedy, 1961), and (2) trans- 
fer of phosphocholine head group from phosphatidylcholine 
(PC) to PC-derived DAG by phosphatidylcholine:diacylglycerol 
cholinephosphotransferase (PDCT) (Lu et al., 2009). Another 
pathway that involves the conversion of monoacylglycerol (MAG) 
to DAG by monoacylglycerol acyltransferase (MGAT) has been 
well characterized in animals (Yen et al, 2002; Cao et al., 2003). 
Although the presence of MAG in plant tissues has been reported 
(Perry and Harwood, 1993; Yang et al., 2002), it is only believed 
to serve as a precursor for cuticular wax biosynthesis (Li-Beisson 
et al, 2013). Recently a plant MGAT gene and an OLEOSIN 
(0LE03) isolated from peanut were shown to enhance TAG 
accumulation when over-expressed in Saccharomyces cerevisiae 
(Parthibane et al., 2012; Vijayaraj et al, 2012). In addition, the 
concept of oil increase in vegetative tissue via a MAG pathway 
was recently demonstrated by (Petrie et al., 2012). The authors 
showed that heterologous expression of mouse MGATs {MGATl 
and 2) in Nicotiana benthamiana leaves can lead to utilization 
of MAG for DAG synthesis, thereby significantly increasing TAG 
accumulation. TAG levels were increased 9.2-fold by MGATl 
and 7.3-fold by MGAT2 at 5 days post-infiltration (dpi) in the 
infiltrated leaves. 

In addition to increased TAG, we found that MGATl infiltrated 
leaf areas of N. benthamiana turn yellow at 3 dpi and exhibit 
senescence-like symptoms by 5 dpi. These symptoms are specific 
to MGATl infiltration and can be reduced by MAG supplementa- 
tion or inhibition of MGAT activity. In order to better understand 
the molecular mechanisms underlying the observed symptoms 
and TAG increase, we undertook a global gene expression study 
of MGATJ -mediated responses in the N. benthamiana transient 
leaf infiltration system. This system has been proven particularly 
useful for the metabolic engineering of lipid biosynthesis path- 
ways (Wood et al., 2009; Petrie et al., 2010; Vanhercke et al., 
2013b) and for functional confirmation of novel genes (Dussert 
et al., 2013). Availability of genetic resources wUl further broaden 
the scope for this model plant to understand the consequences 
of gene/pathway manipulation. Efforts in this direction have 
been initiated with the draft genome being recently published 
(Bombarely et al., 2012) and the transcriptome sequenced from a 
pooled library representing nine different tissues (Nakasugi et al., 
2013). Here we describe the lipid pathway related transcripts of 
N. benthamiana leaf tissue and employed RNA-seq to identify 
genes differentially expressed in MGATl infiltrated leaves. This 
study provides a comprehensive overview of cellular metabolic 
processes caused by the utilization of MAG for TAG synthesis 
and associated senescence-like response resulting from transient 
expression of mouse MGATl in N. benthamiana. 

RESULTS AND DISCUSSION 

PHYSIOLOGICAL AND BIOCHEMICAL RESPONSES OF MGATl 
INFILTRATION 

Transient expression ofM. musculus MGATs (MGATl or MGAT2) 
is known to increase TAG accumulation by 5 days post-infiltration 



(dpi) in N. benthamiana leaves (Petrie et al., 2012). However, the 
MGATJ -infiltrated areas showed cell-death by 5 dpi (Figure lA), 
with senescence-like symptoms being visible from 3 dpi onwards 
(Figure lA). Surprisingly, the MGATl infiltrated leaves show no 
senescence symptoms even at 5 dpi. Differences in expression pat- 
terns, functions and biochemical activities between mammalian 
MGATs have been reported (Cao et al., 2003, 2007). Further, we 
have shown previously that MGATl displays considerable DGAT 
side-activity with 58% of the labeled DAG converted to TAG by 
MGATl in contrast with only 6% conversion by MGAT2 (Petrie 
et al., 2012). Hence, we focused this study on understanding 
the relationship between transient MGATl expression and the 
observed senescence-like responses. As a first step, we measured 
the TAG levels in MGATl infiltrated leaves before the onset of 
senescence symptoms. Lipid analysis on MGATl infiltrated leaves 
at an early time point of 2 dpi showed a 4-fold increase in TAG 
levels (Figure IB). In addition, DAG levels were also increased 
while a clear decrease in MAG levels was observed indicating its 
conversion to DAG. Since TAG accumulation occurs at 2 dpi and 
senescence symptoms do not start to appear, this early time point 
was deemed suitable for gene expression studies. 

DENOVO TRANSCRIPTOME ASSEMBLY AND ANNOTATION 

RNA was extracted from the leaves of 6-week old un-infiltrated 
plants for de novo transcriptome assembly. lUumina sequencing 
followed by processing of raw reads produced a total of 55 mil- 
lion clean reads (Table lA). Assembling of the reads using Trinity 
(Grabherr et al., 2011) produced 92,857 non-redundant contigs 
out of which 77,679 were singletons and 15,178 were distinct clus- 
ters (Table lA). The length distribution of contigs is shown in 
Supplementary Figure 1. The mean contig length was 503 bp with 
a N50 weightage of 706 for the assembly (Table lA). The con- 
tig sequences can be downloaded from CSIRO data access portal 
(http://dx.doi.org/ 1 0.4225/08/524ABB9C42D9 1 ) . 

The 92,857 contigs were searched against various pub- 
lic databases using BLASTX (le^^). The number of tran- 
scripts annotated by each database is summarized in Table IB. 
The NCBI non-redundant (nr) database (ftp://ftp.ncbi.nlm.nih. 
gov/blast/db/) search retrieved annotations for 48,952 contigs. 
Function annotations for 31,572 contigs were obtained from 
Swiss-Prot (Magrane and Consortium, 2011) and "Clusters of 
Orthologous Groups of Proteins" (COG; http://www.ncbi.nlm. 
nih.gov/COG) database search annotated 13,972 contigs. Gene 
Ontology (GO) annotation terms were retrieved for 23,568 con- 
tigs using Blast2GO (Conesa et al, 2005). Thus, a total of 49,180 
contigs were annotated from all the databases indicating that 
nearly 52% of the transcripts were annotated. The list of con- 
tigs and their annotations from each database are described 
in Supplementary data 1. The COG functional class distribu- 
tion and categorization into the GO categories are illustrated 
in Supplementary Figure 2. For metabolic pathway annotations, 
the contigs were searched against Kyoto Encyclopedia of Genes 
and Genomes database (KEGG; http://www.genome.jp/kegg/). 
The enriched pathways and number of transcripts related to 
each pathway are illustrated in Supplementary Figure 3. About 
34.5% (32,030) of transcripts were annotated to the KEGG 
pathways, the majority of which were under the broad cate- 
gory of metabolic pathway (21.28%). Biosynthesis of secondary 
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FIGURE 1 I Physiological and biochemical responses of 
monoacylglycerol acyltransferase {MGATI) transient expression in 
Nicotiana benthamiana leaves. (A) Senescence symptoms observed in 
leaf tissue infiltrated witfi MGATI at 3 days post-infiltration (dpi) and 5 dpi. 
Infiltrations were performed along with p19 viral suppressor protein as 
control (Petrie et al., 2012). (B) Lipid analysis showing monoacylglycerol 
(MAG), diacylglycerol (DAG) and triacylglycerol (TAG) levels in p19 and 
p19+MG/477 infiltrated leaf spots, harvested at 2dpi. Data shown are 
average of three replicates and error bars represent standard error. 



metabolites (9.62%), plant hormone signal transduction (6.19%) 
and plant-pathogen interaction (5.68%) were the other major 
pathway categories (Supplementary Figure 3). From the path- 
way analysis, N. benthamiana contigs related to lipid pathways are 
identified and summarized in Table 2. 

The identified transcripts and their annotations were 
confirmed by randomly BLASTing a number of contig 
sequences against the N. benthamiana transcriptome database 
(http://sydney.edu.au/science/molecular_bioscience/sites/bentha- 
miana/). No disagreement in annotation was observed between 
the matched sequences from both studies. 



Table 1 | Summary of de novo assembling and annotation of 
N. benthamiana contigs. 



SUMMARY OF N. BENTHAMIANA TRANSCRIPTOME SEQUENCES 
GENERATED BY ILLUMINA HiSeqTM2000 

Total raw reads 58,558,070 
Total clearn reads 55,100,052 
Total contigs 92,857 
Singletons 77679 
Distinct clusters 15,178 
Mean length 503 
N50 706 



B 



Database 



No. of N. benthamiana contigs annotated 



Non-redundant (nr) 

SwissProt 

KEGG 

COG 

GO 

Total contigs annotated 



48,952 
31,572 
32,030 
13,972 
23,568 

49,180 



(A) Summary of sequencing and de novo assembling output. (B) Summary of N. 
benthamiana contigs annotated by different databases. 

Table 2 | Summary of N. benthamiana contigs related to different 
lipid pathways identified by KEGG pathway analysis. 



KEGG Pathway 



Number Percentage Pathway ID 



81 


0.33 


koOOOOl 


121 


0.50 


ko00071 


543 


2.24 


ko00564 


345 


1.43 


ko00565 


115 


0.48 


koooeoo 


63 


0.26 


ko00604 


114 


0.47 


ko01040 



Fatty acid biosynthesis 
Fatty acid metabolism 
Glycerophospholipid metabolism 
Ether lipid metabolism 
Sphingolipid metabolism 
Glycosphingolipid 
biosynthesis — ganglio series 
Biosynthesis of unsaturated fatty 
acids 



IDENTIFICATION OF LEAF ACYL-LIPID METABOLISM GENES OF 
N. BENTHAMIANA 

Acyl-lipids include a wide range of lipid classes with diverse func- 
tions in cell membranes and sub-cellular organelles (Li-Beisson 
et al., 2010). Identification of acyl-lipid genes of N. benthami- 
ana is important in the context that this model plant can be 
used extensively for manipulating lipid pathways in leaf tissue. 
The N. benthamiana contig sequences were mapped against the 
Arabidopsis TAIRIO protein sequences (http://www.arabidopsis. 
org/) using BLASTX (le^^) in CLC Genomics Workbench 5 
(http://www.clcbio.com). Identified contigs with an Arabidopsis 
thaliana homolog were analyzed against acyl-lipid metabolism 
gene list obtained from "The Arabidopsis Acyl-Lipid Metabolism" 
website (ARALIP; http://aralip.plantbiology.msu.edu). A total of 
1648 acyl-lipid contigs of N. benthamiana leaf were identified. 
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FIGURE 2 I Identification of leaf acyl-lipid metabolism transcripts of 
N. benthamiana. Contigs mapped to A. thaliana sequences were searched 
against "The Arabidopsis Acyl-Lipid Metabolism" website (ARALIP; 
http://aralip.plantbiology.msu.edu) to identify acyl-lipid related transcripts. 
Distribution of identified transcripts into different acyl-lipid metabolism 
categories and number of contigs in each category is illustrated. 



The distribution of these transcripts into 23 different lipid- 
related function categories and transcript count for each cate- 
gory are shown in Figure 2. In plants, the majority of the seed 
oil is stored as TAG. The N. benthamiana transcripts that are 
involved in TAG biosynthesis were identified and are described 
in Supplementary Table 1. A list of all the contigs identified 
with a role in acyl-lipid metabolism is shown in Supplementary 
data 2. 

TRANSCRIPTIONAL RESPONSES OF MGAT1 INFILTRATION 

In order to understand the molecular mechanisms underly- 
ing MGATi -mediated TAG accumulation and senescence-like 
responses, RNA-seq analysis was performed on MGATI infil- 
trated vs. control infiltrated leaves. The early time point of 2 dpi 
showing TAG accumulation but not senescence symptoms along 
with an additional time point of 1 dpi was selected for iden- 
tifying differentially expressed genes. RNA was extracted from 
the control (pl9) and MGATJ -infiltrated leaf areas, pooled from 
at least 10 different leaves for each treatment at each time 
point, and lUumina sequencing reads were obtained. The expres- 
sion levels were estimated as Reads Per Kilobase of Exon per 
Million mapped reads (RPKM) by mapping against the de novo 
assembled transcriptome (92,857 contigs) in CLC Genomics. 
Differentially expressed transcripts were identified according to 



the RNA-seq workflow in CLC Genomics and also by using a R 
package DESeq (Anders and Huber, 2010). No significant dif- 
ferences were observed in the fold changes derived from both 
tools. Gene lists derived from DESeq containing p-values (cut-off 
p < 0.05) adjusted for multiple testing by Benjamini-Hochberg 
procedure were used for further analysis. After applying a fold 
change (EC) cut-off of >2/<— 2 (p < 0.05), 178 genes were up- 
regulated and 134 genes were down-regulated at 1 dpi. At 2 dpi, 
846 genes were up-regulated and 2434 genes down-regulated 
(Supplementary Eigure 4A). A. thaliana homologs were identi- 
fied for 152 (1 dpi) and 482 (2 dpi) up-regulated genes, and 103 
(1 dpi) and 1890 (2 dpi) down-regulated genes. It is note-worthy 
that, similar to MGATI response at 2 dpi, several senescence 
studies noted a larger number of genes being down-regulated 
than up-regulated genes under senescing conditions (Troncoso- 
Ponce et al., 2013). A scatter plot using the log2 EC values 
of 4 genes estimated by RNA-seq and qRT-PCR is shown in 
Supplementary Eigure 4B. Eurther, the expression of MGATI 
in infiltrated leaves was confirmed by RT-PCR (Supplementary 
Figure 4C). 

GENE ONTOLOGY (GO) AND PATHWAY (MapMan) ANNOTATIONS OF 
DIFFERENTIALLY EXPRESSED GENES 

Further analysis was performed on contigs having an A. thaliana 
homolog by taking advantage of the available resources. Since 
major changes in gene expression were observed at 2 dpi only, this 
category of genes was selected for further analysis. After remov- 
ing duplications in annotations, 1666 transcripts were found to 
be changed at 2 dpi, out of which 1333 were down-regulated 
and 333 were up-regulated. The up- and down-regulated genes 
of 2 dpi were analyzed for enriched GO terms {p < 0.05) using 
Classification SuperViewer tool at The Bio-Analytic Resource 
(BAR) for Plant Biology (http://bar.utoronto.ca), as illustrated in 
Figure 3. The up-regulated genes were highly enriched in stress 
response, abiotic or biotic stimulus and signal transduction bio- 
logical responses (Figure 3A). This indicates that, in addition to 
accumulation of increased TAG, MGATI infiltrated areas were 
also experiencing and responding to stress. A significant num- 
ber of down-regulated genes were located in plastid, chloroplast 
and ceU-waU (Figure 3B), suggesting that cellular processes local- 
ized to these organelles could be negatively impacted by MGATI . 
Indeed the electron transport and energy pathways were highly 
enriched in down-regulated genes indicating a suppressive effect 
on these processes by MGATI . 

METABOLISM OVERVIEW 

To provide an overview of MGATI effect on individual metabolic 
pathways, significantly changed genes (>2/<— 2-fold; p < 0.05) 
along with their fold changes were mapped on pathway bins 
of MapMan 3.5. 1R2 (http://mapman.gabipd.org). Metabolism 
overview revealed that the majority of MGATI responsive genes 
map to cell- wall, lipids and PS bins (Figure 4A). Most of 
the cell-wall related genes including cellulose synthases (CESA 
2-6), cellulose synthase like (CSL) genes, and arabinogalac- 
tan proteins (AGPs) such as salt overly sensitive 5 {SOS 5) 
involved in cell-wall biosynthesis (Xu et al., 2008) were down- 
regulated. CeU-waU modifying genes like expansins, xyloglucan 
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FIGURE 3 I Gene ontology (GO) annotations of differentially 
expressed genes (DEG). Classification SuperViewer tool at The 
Bio-Analytic Resource (BAR) for Plant Biology (http;//bar.utoronto.ca) 



was used. (A,B) GO term annotation of up-regulated (A) and 
down-regulated genes (B). Significantly enrictned terms are 
highlighted in bold. 



endotransglycosylase (XETs) and pectin esterases were also down- 
regulated. Down-regulation of cell-wall synthesis and modifica- 
tion genes indicates a possible disruption in cell-wall integrity. 
Genes involved in carbohydrate (CHO) metabolism were also 
affected by MG ATI (Figure 4A). While many major CHO genes 
related to sucrose and starch metabolism were down-regulated, a 
cell-wall invertase (AtcwINV2; At3g52600) and a vacuolar inver- 
tase ('Atlg62660) were up-regulated. Decline of photosynthesis 
(PS) initiates senescence and a role for invertases in this process 
has been demonstrated (Dickinson et al, 1991; Ding et al., 1993). 
Minor CHO related genes including trehalose-6-phosphate phos- 
phatases (TPPs) and two trehalose-phosphate synthases (TPS 
5 and 11) were also up-regulated by MGATl (Figure 4A). An 
overview of MGATl effect on PS is illustrated in Figure 4B and 
corresponding genes are summarized in Supplementary Table 
2. In agreement with the GO analysis showing down-regulated 
terms of electron transport, energy and plastid, PS showed largest 
down-regulation effect by MGATl, as compared to other path- 
ways. The collapse of PS, cell-wall, major CHO metabolism and 
electron transport/ ATP synthesis are major indicators of senes- 
cence induction that can cause the visible symptoms observed 
at 3 dpi. 



GENES INVOLVED IN MGATl INDUCED TAG BIOSYNTHESIS 

The N. benthamiana acyl-lipid transcripts identified in the study 
were analyzed for differentially regulated genes and lipid processes 
induced by MGATl infiltration. Expression of 94 transcripts 
mapped to ARALIP genes were altered by MGATl expression, 
of which only 12 were up-regulated (Supplementary Table 3). 
The only lipid process enriched with up-regulated genes was 
TAG biosynthesis, while genes involved in FA synthesis and 
modification, TAG degradation, FA elongation, wax biosynthe- 
sis, cutin and suberin synthesis and transport were all found to 
be down-regulated (Figures). These responses are indicative of 
the induction of TAG biosynthesis machinery by MGATl expres- 
sion while trancripts involved in FA degradation, synthesis and 
transport of wax, suberin and cuticular structures appear to be 
suppressed. Two TAG biosynthesis related genes, diacylglycerol 
acyltransferase (DGATl ) and phosphatidylcholine:diacylglycerol 
cholinephosphotransferase (PDCT) were up-regulated. FA degra- 
dation genes down-regulated included two MAG lipases (MAGL) 
and a TAG lipase (TAGL). PDCT and the transiently expressed 
MGATl are involved in synthesizing DAG from PC and MAG, 
respectively. DAG is the only substrate of DGATl for conver- 
sion to TAG. Transient expression of MGATl and the subsequent 
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FIGURE 4 I MapMan pathway overview analysis of DEC (A) Overview 
of MGAT1 effect on transcripts related to different pathways generated by 
"Metabolism overview" analysis in MapMan. Fold changes (Log2) are 
indicated by the colour scale. Distribution of transcripts into different 
pathways and their expression trend in each pathway is illustrated. 



(B) Effect of MGATl on photosynthesis (PS) related transcripts. All PS 
related genes were down-regulated. The number of transcripts, their 
distribution into sub-classes of PS and expression trend are illustrated. The 
17 sub-classes and the genes in each sub-class are described in 
Supplementary Table 2. 
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up-regulation of PDCT and DGATl could explain the increased 
DAG and TAG levels detected in MGATl infiltrated areas in com- 
bination with depletion in MAG levels compared to the negative 
control (pl9). Down-regulation of MAG and TAG lipases could 
result in reduced neutral lipid turnover thereby contributing to a 
net increase in TAG levels. 

EXPRESSION PROFILES OF MGAT1 RESPONSIVE LIPID-RELATED 
GENES DURING SENESCENCE 

Recently, Troncoso-Ponce et al. (2013) have developed a database 
by compiling and plotting the senescence responsive profiles of 
lipid-related genes from three independent studies (Schmid et al, 
2005; Van Der Graaff et al, 2006; Breeze et al, 2011). We chose 
Breeze et al. (2011) and Van Der Graaff et al. (2006) studies only 
to compare the MGATl responsive profiles of the lipid-related 
genes as Schmid et al. study was not specifically directed at senes- 
cence. Breeze et al., analyzed expression pattern of A. thaliana 
genes from 19 days after sowing (DAS) to 39 DAS. Van der Graaff 
et al. studied expression profiles of developmental leaf senes- 
cence (NS), darkening-induced senescence in intact leaves (DIS) 
and senescence in dark-incubated detached leaves (DET) at dif- 
ferent stages. Out of the 12 lipid-related genes identified to be 
up-regulated by MGATl, 11 genes were also up-regulated dur- 
ing senescence in both studies (Figure 6A). A large number (69) 
of lipid-related genes were down-regulated by MGATl. While 
no clear trend emerged, the majority of these genes were also 
down-regulated in both Breeze et al. (78%) and Van der Graaff 
et al. (52%) senescence studies (Supplementaty Figure 5). These 
results indicate that transcriptional changes caused by MGATl 
infiltration are similar to those induced by senescence as discussed 
below. 

One of the characteristic responses of lipid-related genes dur- 
ing senescence is a co-ordinated decrease in the transcripts for 
FA synthesis and desaturation, while the glycerolipid assembly 
transcripts including DGATl increase in expression (Troncoso- 
Ponce et al., 2013). In line with this, transcripts involved in FA 
synthesis such as KASI, KASII, P-ketoacyl-CoA synthase (KCS) 
and long-chain acyl-CoA synthase {LACS2), as well as FAD3 
and 5 involved in FA desaturation were down-regulated by 
MGATl, while DGATl and PDCT involved in glycerolipid assem- 
bly were up-regulated. In contrast to the up-regulation of some 
glycerol-3-phosphate acyltransferases {GPAT9 and GPATl) dur- 
ing senescence (Troncoso-Ponce et al., 2013), GPAT4 required 
for MAG accumulation was down-regulated by MGATl. This 
finding indicates that in addition to utilizing MAG for conver- 
sion to DAG, MGATl also affects MAG synthesis by decreasing 
GPAT4 expression which can result in further reduction of MAG 
levels in MGATl infiltrated leaves. Up-regulation of genes encod- 
ing various lipases except for phospholipase D (PLD) and genes 
involved in ^-oxidation pathway is a feature of lipid turnover 
during natural senescence (Thompson et al, 1998; Troncoso- 
Ponce et al., 2013). Somewhat unexpectedly, transcripts of two 
putative MAG lipases and a TAG lipase were down-regulated 
by MGATl expression, while no change was observed for P- 
oxidation genes. In addition, PLD-P and PLD-S transcripts 
were up-regulated by MGATl. Several putative ABC transporters 
involved in FA transport were also down- regulated by MGATl. 



Blocking of FA transport to peroxisomes by mutating an ABC 
transporter resulted in increased TAG levels during senescence 
(Slocombe et al, 2009). Thus, whUe several MGATl induced 
lipid-responses are similar to senescence responses including 
increase in DGATl transcripts and higher TAG levels, it appears 
that MGATl also blocks FA turnover resulting in reduced TAG 
turnover. 

COMPARATIVE ANALYSIS OF MGAT1 AND SENESCENCE RESPONSIVE 
GENE EXPRESSION 

To understand the relation between MGATl expression and 
senescence, the entire set of MGATi -responsive genes was com- 
pared with senescence response gene lists from both Breeze et al. 
(6323 genes) and Van der Graaff et al. (4109 genes) studies. The 
expression of 466 out of 1666 (27.9%) MGATi -responsive genes 
were also changed in both studies. While 314 genes showed a 
specific overlap with Breeze et al., 191 MGATi -responsive genes 
overlapped specifically with Van der Graaff et al. As expected, 
a larger proportion (2532) of genes was common between the 
two senescence studies. About 47.1% (695) of MGATi -responsive 
genes were specific to MGATl and were unaffected in either study 
(Figure 6B). Interestingly, this subset includes transcription fac- 
tors MYB5 and TRANSPARENT TESTA GLABRA (TTG2), found 
to be down-regulated by MGATl. In A. thaliana, MYB5 and 
rrG2 are key components in regulation of mucilage biosynthe- 
sis (Western et al., 2004). While myb5-l mutants exhibit loss of 
mucilage and enhanced oil accumulation (Shi et al, 2012), rrG2 
is required for the expression of GLABRA2 (GL2), a negative reg- 
ulator of seed oil accumulation (Shen et al., 2006; Shi et al, 2012). 
Taken together, it is apparent from these results that in addition 
to a good overlap with senescence related genes, MGATi -specific 
changes were also abundant. 

In order to perform comparisons based on expression pattern 
and associated functions, cluster analysis was performed on 
gene lists from all three studies using MapMan. Genes from 
each study were categorized into up- and down-regulated 
and given a value of "2" and "—2" respectively. Genes that 
did not change in a given condition or gene list were given 
a value of "1." Using eucledian distance metric, 15 clusters 
were generated (Supplementary Figure 6). Interestingly, 216 
genes were down-regulated in all studies and grouped to 
Cluster Al. Out of these 24% belong to PS. Genes related to 
ceU-wall, lipid and protein metabolism, secondary metabolism 
and RNA regulation were also present in this cluster. Clusters 
A2-5 consists of genes down-regulated by MGATl and in one 
or more other studies. Clusters A3 (56 genes) and A4 (30 
genes) together represent the second largest set of genes. This 
set includes many transcription factors and signaling kinases 
along with FA synthesis genes BiASII and LACS (A3). These 
profiles point to the similarities of MGATl responses to nat- 
ural senescence in affecting PS and lipid pathways. Cluster Bl 
genes were up-regulated in almost every condition. This subset 
includes SRG3, PR4 (PATHOGENESIS-RELATED 4), heat shock 
protein {HSP17.8-C1), HSP70, SIP4 iSOS3-INTERACTING 
PROTEIN 4) and CML37 (CALMODULIN LIKE 37) with known 
roles in biotic and abiotic stresses. The key gene for ethylene 
biosynthesis, ACS2 (Atlg01480) is also present in this cluster. 
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FIGURE 5 1 Identification of MGAT1 DEGs involved in acyl-lipid 
metabolism. Distribution of genes differentially expressed by MGATl into 
different categories of acyl-lipid metabolism is illustrated. Down-regulation 


(blue) of transcripts in majority of classes including FA elongation, cutin, 
suberin and wax biosynthesis can be seen, while TAG biosynthesis contains 
more up-regulated (green) genes. 



Up-regulation of these genes in all conditions is indicative of 
stress and senescence related responses and increased synthe- 
sis of the senescence hormone ethylene. Genes up-regulated 
by MGATl but either down-regulated or unchanged in one 
or more of the other conditions were grouped in clusters 
B2 to B5. This set includes a senescence marker gene YLS9 
{YELLOW-LEAF-SPECIFIC GENE 9) that is responsive to late 
senescence, dark, ethylene and ABA (Yoshida et al., 2001), 
PLDBETAl {PHOSPHOLIPASE D BETA 1) [cluster B3]; the ABA 
biosynthesis gene, NCED3 {NINE-CIS-EPOXYCAROTENOID 
DIOXYGENASE 3), two Ethylene signaling genes EBF2 



{EIN3-BINDING F BOX PROTEIN 2) and RRTFl (REDOX 
RESPONSIVE TRANSCRIPTION FACTOR 1), and the tran- 
scription factor STZ (salt tolerance to zinc finger) [Cluster B4]. 
STZ is a key factor controlling the expression of senescence 
related NAC transcription factors (Breeze et al., 2011). Cluster 
CI consists of 5 genes specifically up-regulated by MGATl which 
include a cell-wall synthesis related gene MATERNAL EFFECT 
EMBRYO ARREST (MEE31) and a peroxidise (PER12). Genes 
down-regulated specifically by MGATl were grouped in clusters 
C2 (31 genes) and C3 (30 genes). Among the C3 genes, 7 were 
transcription factors which include two ethylene signaling genes 
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FIGURE 6 I MGATl and senescence responsive transcripts. 

(A) Expression profiles of MGATl responsive lipid related genes during 
senescence. Senescence profiles from two earlier studies (Van Der 
Graaff et al., 2006; Breeze et al., 2011) were downloaded from 
Troncoso-Ponce et al. (2013) database and shown. Eleven out of twelve 



lipid-related genes up-regulated by MGATl were also up-regulated during 
progression of senescence shown as days after sowing (DAS). 
(B) Comparative analysis of MGATl and senescence responsive 
transcriptome. Venn diagram showing gene overlap between MGATl, 
Breeze et al. and Vander Graff et al. 



ETHYLENE-RESPONSIVE ELEMENT BINDING PROTEIN 
{ATEB, At3gl6770) and ETHYLENE-INSENSITIVE 3 (EIN3). 
Genes down-regulated by MGATl and in DET were grouped in 
cluster C4. Important in this cluster was an ABA biosynthesis 
gene, NCED4 that was down-regulated in NS as well. These 
results point to senescence-like gene expression changes in 
response to MGATl infiltration, the majority of which affect PS. 
Similarities in phytohormone-related gene expression patterns 
were also apparent. 



MGATl AND LEAF SENESCENCE 

To confirm the role of MGATl in the observed senescence, 
N. benthamiana leaves were infiltrated with four different 
concentrations of Agrobacterium tumefaciens (AGLl) containing 
MGATl. A clear correlation between MGATl concentration and 
leaf yellowing was observed at 3 dpi (Figures 7Ac-f). No signifi- 
cant yellowing was visible at 0.03 OD concentration (Figure 7Ac), 
while yellowing appeared at 0.0625 OD (Figure 7Ad) and was 
much more pronounced at 0.125 (Figure 7Ae) and 0.25 OD 
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(Figure 7 AT). Infiltrations with A. tumefaciens cultures contain- 
ing pl9 at OD 0.125 and 0.375 did not result in a senescence 
phenotypic response (Figures 7Aa,b). Although the above results 
indicated a clear role of MGATl in the induction of senescence- 
like responses, we further confirmed this hypothesis by blocking 
the MGATl activity. The zwitterionic detergent, CHAPS has been 
used extensively to inhibit the MGAT and DGAT activities of 
MGAT enzymes (Cao et al, 2007; Vijayaraj et al, 2012). Co- 
infiltration of 1 mM CHAPS along with MGATl did not show 
senescence symptoms in the infiltrated area (Figures 7Bg) while 
symptoms were observed in the MGATl infiltrated areas of the 
same leaves (Figure 7Bd). These results clearly show that the 
activity of infiltrated MGATl is required for observed senescence- 
like symptoms. 

ROLE OF MAG IN MGAT1-INDUCED SENESCENCE-LIKE RESPONSE 

Although we confirmed the role of MGATl in the induction of 
senescence-like symptoms, the involvement of MAG depletion 
or DAG/TAG accumulation was still not clear. Since accumu- 
lation of DAG or TAG has not always been associated with 
senescence, we hypothesized that decrease in MAG levels by 
MGATl activity could be involved. However, since it is known 
that co-infiltration of DGATl together with MGAT2 leads to 
the efficient conversion of accumulated DAG to TAG by the co- 
expressed DGATl (Petrie et al., 2012), we first tested the effect 
of DGATl co-infiltration with MGATl on observed senescence. 
As shown in Figure 7B, no significant changes in visible senes- 
cence were observed between MGATl (Figure 7Bd) infiltration 
and co-infiltration with DGATl at two different concentrations 
of 0.125 and 0.25 OD (Figures 7Be-f). DGATl itself did not 
cause any senescence (Figure 7Bc). These results, although not 
conclusive, indicate that DGATl activity had no effect on senes- 
cence. The fact that the observed senescence was comparable 
between MGATl and MGATl+DGATl pointed to a relation 
between reduced MAG levels and senescence. Indeed, supple- 
mentation with an 18:1 sn-2 MAG (2-Oleoylglycerol) at dif- 
ferent concentrations (10-60 |iM) during infiltration alongside 
MGATl rescued the observed senescence symptoms in a con- 
centration dependent manner (Figures 7Cc-f) as compared to 
MGATl only (Figure 7Cb). To further confirm the role of MAG 
depletion in senescence, GPAT4 gene encoding a sn-2 specific 
acyltransferase involved in MAG synthesis was co-infiltrated with 
MGATl. Similar to MAG supplementation, GPAT4 co-infiltration 
(Figure 7Dd) reduced the severity of senescence symptoms 
induced by MGATl (Figure 7Db). A tight correlation between 
observed senescence and supplementation by both exogenous 
and genetic modes are indicative of a role for MAG depletion in 
MGATl induced senescence. It is noteworthy that GPAT4 was also 
down-regulated by MGATl. 

ION-CONDUCTIVITY MEASUREMENT OF MGATl INDUCED 
SENESCENCE-LIKE RESPONSE 

In order to quantify the cell damage, electrical conductivity 
of infiltrated leaf areas was measured as described in meth- 
ods. Leaves were infihrated with MGATl, MGATl + 60|xM 
MAG, MGATl + GPAT4 and MGATl + 1 mM CHAPS alongside 
appropriate controls for each treatment. Leaf discs harvested 



from infiltrated leaves were floated overnight in sterile water 
for measuring conductivity. Similar to the visual observations, 
MGATl infiltrated leaves showed highest conductivity (30|xS). 
Both MAG (60 [lM) supplementation and GPAT4 co- infiltration 
reduced the conductivity by 24.5% and 29.7% respectively, 
whereas a 34% reduction was observed in presence of 1 mM 
CHAPS (Figure 7E). All controls including pl9 showed a mean 
conductivity of 15 (xS. These results show that MAG depletion 
due to conversion to DAG by MGATl can induce senescence- 
like symptoms that can be rescued by providing additional MAG 
substrate. Measurement of TAG, DAG and MAG levels under 
different senescence rescue treatments should provide further 
biochemical evidence for the role of MAG in senescence. 

CONCLUSIONS 

We have shown that recruiting MAG for TAG synthesis via 
MGATl transient expression results in increased TAG accumula- 
tion as early as 2 dpi and also induces senescence-like response 
in N. benthamiana infiltrated leaves. A de novo transcriptome 
data set was generated from leaf tissue to study the underlying 
molecular changes. This data set allowed identification of genes 
involved in lipid pathways of the model plant N. benthamiana 
and can serve as a resource for future transcriptomic studies 
in this system. Using this transcriptome as reference we have 
identified genes differentially expressed by MGATl infiltration. 
Similarities in transcriptional responses to MGATl expression 
and development of senescence were identified by comparative 
analysis. These similarities include a majority of down-regulated 
genes, down-regulation of entire set of photosynthesis related 
genes, activation of senescence hormone responses, decrease in 
FA synthesis and desaturation related transcripts. On the other 
hand, genes involved in glycerolipid assembly such as PDCT 
and DGATl were found to be up-regulated. These results pro- 
vide a mechanism for MGAT-mediated TAG synthesis in which 
transient expression of MGATl and subsequent up-regulation 
of PDCT results in enhanced DAG synthesis from MAG and 
PC, respectively. Subsequently, DAG is efficiently converted to 
TAG due to (a) the co-ordinated up-regulation of DGATl and 
(b) the DGAT side-activity of MGATl (Figures). In addition, 
interesting MGATl responses were identified including the down- 
regulation of transcription factors MYB5 and TTG2 involved in 
the regulation of mucilage biosynthesis. The senescence related 
responses of MGATl include dramatic decrease in PS-related 
genes, down-regulation of cutin, suberin and wax biosynthe- 
sis genes including GPAT4 (Figure 8). Further, we demonstrate 
that MGATl inhibition, addition of MAG or co- infiltration of 
GPAT4 can reduce the observed senescence as indicated by a 
decrease in ion-conductivity. Taken together, these findings sug- 
gest that MAG conversion to DAG by MGATl leads to MAG 
depletion affecting diverse biological processes via gene expres- 
sion changes thereby initiating a senescence-like response. DAG 
has been suggested as a signaling molecule (Dong et al., 2012). 
While its role in MGATl mediated senescence-like response 
cannot be ruled out, the N. benthamiana system will allow 
further characterization of the role of different DAG molecu- 
lar species under senescence-promoting and -rescue conditions. 
Finally, MGATl transient expression could provide a valuable 
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FIGURE 7 I Effect of MGATl expression and IVIAG depletion on observed 
senescence symptoms. (A) Leaf spots infiltrated with two different 
concentrations of pl9 lOD 0.125 (Aa) and OD 0.375 (Ab)l and four different 
concentrations of pl9 + MGATl 10.03 OD (Ac), 0.0625 OD (Ad), 0.125 OD 
(Ae) and 0.25 OD (Af)]. The severity of visible symptoms is correlated to 
MGATl concentration. (B) Co-infiltration of CHAPS or diacylglycerol 
acyltransferase {DGATl) with MGATl. No senescence symptoms are 
detected in the control treatments— pi 9 (a), pl9 + 1 mM CHAPS (b) and 
pl9 + DGATl (c). Senescence induced by MGATl (d) is abolished by 
including 1 mM CHAPS in infiltration mix containing MGATl (g). No effect 
was observed when DGATl was co-infiltrated [0.125 OD (e) and 0.25 OD (f)]. 
(C) MAG supplementation: 2-Oleoylglycerol was added to the infiltration 



buffer to obtain IOm-M (c), 20m.M (d), 40|iM (e), or 60|iM (f) final 
concentrations. Senescence was considerably reduced in the presence of 
60|iM MAG (f) in comparison to MGATl alone (b). Increase in MAG 
concentration was accompanied by a decrease in visible senescence (c-f). 
(D) Co-infiltration of a glycerol-3-phosphate acyltransferase (AtGPAT4) with 
MGATl (d) considerably reduced visible senescence as compared to the 
MGATl control (b). Infiltration of AtGPATA by itself (c) did not induce 
senescence. (E) Ion-conductivity measurement of MGATl induced 
senescence. Ion conductivity of infiltrated leaf discs was measured by 
floating in water Water with no leaf discs floated (H2O) and un-infiltrated (WT) 
leaf discs served as controls. Data shown are average of three replicates and 
error bars represent standard error 
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FIGURE 8 I Schematic overview of MG4r7-induced TAG increase and 
senescence-lil<e response. Up- (red) and down-regulated (hue) genes 
are highlighted. Transient expression of MGATl and subsequent 
up-regulation of PDCT possibly result in enhanced DAG synthesis from 
MAG and PC, respectively. DAG is subsequently converted to TAG due 
to the co-ordinated up-regulation of DGATl and DGAT side-activity of 
MGATl. Down-regulation of TAG and MAG lipases {TAGL and MAGL) 



indicate decreased FA degradation and increased TAG retention. 
Down-regulation of cutin, wax and suberin biosynthesis transcripts 
including GPAT4 and diversion of MAG for TAG synthesis is likely to 
affect biosynthesis of surface lipids. In addition, down-regulation of 
photosynthesis genes and up-regulation of senescence hormone related 
genes indicate a co-ordinated mechanism that might initiate 
senescence-like response. 



tool for the study of surface lipids biosynthesis such as cutin, 
wax and suberin as well as plant-pathogen research due to the 
ability to induce controlled cell-death in the leaves of a model 
plant. 

MATERIALS AND METHODS 

TRANSIENT EXPRESSION IN N. BENTHAMIANA 

All infiltrations of Agrobacterium tumefaciens (AGLl) for tran- 
sient gene expression in N. benthamiana were performed on 
6 week old leaves as described previously (Voinnet et al., 
2003; Wood et al, 2009). RNA extracted from 6 weeks 
old un-infiltrated N. benthamiana leaves was used for de 
novo transcriptome assembly. For differential gene expression 
analysis, MGATl and pl9 genes were transiently expressed. 
Infiltrated areas from at least 10 different leaves were pooled 
after 1 and 2 dpi separately and frozen for RNA extraction. 
10 mM CHAPS (3-[(3-Cholamidopropyl)dimethylammonio]-l- 
propanesulfonate); (Sigma-Aldrich, USA) solution prepared in 
deionized water was included in the infiltration buffer to 
yield a final concentration of 1 mM. For MAG supplemen- 
tation, 2-Oleoylglycerol (Sigma-Aldrich) dissolved in acetoni- 
trUe was added to the infiltration buffer to obtain 10, 20, 



40, or 60 [xM final concentrations. The mouse MGATl and 
A. thaliana DGATl and GPAT4 sequences used in this study 
were described previously (Yang et al, 2002; Petrie et al, 
2012). 

TRANSCRIPTOME ANALYSIS 
lllumina sequencing 

RNA (20 |ig) extracted from frozen leaf tissues using RNeasy 
Plant Mini kit (Qiagen, Hilden, Germany) was used for cDNA 
library preparation and sequencing at Beijing Genome Institute 
(BGI, Shenzhen, China). Briefly, poly (A) mRNA was isolated 
with Oligo(dT) attached beads and subjected to fragmentation. 
First strand cDNA synthesis by random hexamers was followed by 
second-strand synthesis and purification (QiaQuick PGR purifi- 
cation kit, Qiagen), end repair and adapter ligation. Equimolar 
concentrations of single end, 100 bp libraries generated after size 
selection and PGR enrichment were sequenced in a single lane of 
lllumina HiSeq™ 2000 system. 

De novo assembly and annotation 

The RNA-Seq reads generated were initially processed to 
remove the adapter sequences, ambiguous reads (5%) and low 
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quality sequences (Q20). Reads were de novo assembled to 
contigs using Trinity (Grabherr et al., 2011). Singletons and 
clusters with isoforms were identified by building de Brujin 
graphs using Chrysalis module in Trinity. The contigs were 
mapped against public databases including nr (ftp://ftp.ncbi.nlm. 
nih.gov/blast/db/), Swiss-prot (http://www.uniprot.org/), COG 
(http://www.ncbi.nlm.nih.gov/COG) and KEGG (http://www. 
genome.jp/kegg/) using BLASTX (ii-value cutoff le^^) for func- 
tional annotation. Gene Ontology (GO) annotation terms were 
obtained by using Blast2GO (Conesa et al, 2005). The con- 
tigs were aligned against the A. thaliana protein database TAIR 
10 (http://www.arabidopsis.org/) by downloading the sequences 
and BLASTing (cut-off le^^) in CLC Genomics Workbench 5 
(http://www.clcbio.com). The A. thaliana acyl-lipid gene list was 
obtained from "The Arabidopsis Acyl-Lipid Metabolism" web- 
site (ARALIP; http://aralip.plantbiology.msu.edu). The GO term 
enrichment analysis was performed by Classification SuperViewer 
tool at BAR (http://bar.utoronto.ca). 

Differential gene expression analysis 

RNA-Seq reads from the MGATl and pl9 (control) RNA samples 
were loaded on to CLCBio Genomics workbench 5 for expres- 
sion analysis. The expression values of transcripts in each sample 
were measured as RPKM (Mortazavi et al, 2008) in the RNA- 
seq workflow of CLC Genomics. The parameters used include 
maximum number of 2 mismatches and minimum read count 
of 10. Fold change analyses was performed in CLC Genomics 
by setting up a two group comparison experiment. Differential 
expression analysis was also carried using DESeq package (1.14.0) 
in R (2.15.3). The nbinomTest function was used to estimate 
fold changes which employs Benjamini-Hochberg procedure for 
false discovery rate (FDR) control and outputs the fold changes 
and the p-values for statistical significance (Anders and Huber, 
2010). Transcripts with 2-fold up- or down-regulation (p < 0.05) 
and having assigned A. thaliana genes were further analyzed. 
The A. thaliana gene expression data during natural senescence 
(Troncoso-Ponce et al, 2013) was used for comparing the MGATl 
responsive expression profiles at day 2. The MGATl responsive 
gene list (1666 genes). Supplementary data set of Breeze et al. and 
Supplement Table 3 of Vander Graaff et al. were used for compar- 
ative analysis and for cluster analysis in MapMan 3.5. 1R2 (http:// 
mapman.gabipd.org). 

qRT-PCR ANALYSIS 

Total RNA (1 was converted to cDNA using First- Strand 
cDNA synthesis mix (OriGene Inc., Australia). Diluted cDNA 
(0.2X) was used for quantitative real-time PGR by Bio-Rad 
CFX Real-Time System (BIO-RAD, USA) using iQSYBR Green 
Supermrx (BIO-RAD). Primer validation was performed for all 
target genes with five dilutions of template concentration along 
with the N. benthamiana ACTIN and GAPDH control genes. 
Reactions were carried out with initial denaturation at 95°C for 
3min, followed by 35 cycles of 95°C for 10 s, 58°C for 30 s and 
68°C for 30 s. Relative fold differences were determined by A ACj 
method (Livak and Schmittgen, 2001). Primer sequences and 
description of the target and control sequences are provided in 
Supplementary Table 4. 



ELECTROLYTE LEAKAGE 

Six week old N. benthamiana leaves were infiltrated as described 
above. Infiltrated leaf discs (0.5 cm diameter) were harvested 
after 2 days, rinsed twice (2-3 min) with deionized water and 
subsequently floated on 15 mL of deionized water with shaking 
(80 rpm) at room temperature. The electrolyte leakage in the 
solution was measured after 17h of floating using a conductiv- 
ity meter (CD-4301, Lutron). All measurements were done on 
three biological replicates by pooling samples from six plants (two 
plants per replicate). Values are the average of three replicates and 
error bars represent standard error. 

LIPID ANALYSIS 

Leaf discs (2 cm diameter) were cut from 6 different plants (two 
plants per replicate). Total lipids were extracted from freeze dried 
tissue as previously reported (Petrie et al., 2012), run on a TLC 
plate (20 cm^. Silica gel 60, Merck) and developed using a two- 
step solvent system. Samples were first run to 12 cm in chloro- 
form:methanol:acetic acid:water (68:22:6:4, v:v:v:v), followed by 
a second separation in hexane:diethyl ether:acetic acid (70:30:1, 
v:v:v). TAG, DAG and MAG bands were visualized using iodine 
vapor and isolated from the silica gel. Fatty acid methyl esters 
were prepared from each isolated lipid fraction and analyzed by 
GC-FID (7890A GC, Agilent Technologies, Palo Alto, CA, USA) 
equipped with a 30 m BPX70 column (SGE, Austin, TX, USA) 
as described previously (Petrie et al., 2012). Peaks were inte- 
grated with Agilent Technologies ChemStation software (Revion 
B.04.03). 
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